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Abstract 


A nonlinear wave mechanical equation is proposed by inserting an 
imaginary quantum potential into the Schrödinger equation. An explicit 
expression for its solution is given under certain assumptions and it is 
shown that it entails attenuation related effects as non-unitary evolution, 
non-exponential quantum decay and entropy production. In the quantum 
hydrodynamical formulation the existence of circulation effects for the os- 
motic velocity field is established. Finally, a time-invariant equation for 
the probability density is derived, analogous to the tensor Lighthill equa- 
tion in aeroacoustics, which admits both retarded and advanced solutions. 


1 Introduction 


One of the most spectacular aspects of quantum theory is that it can be formu- 
lated in a variety of ways that even appear not to be equivalent at first sight. 
In some occasions the mathematical formulation is identical but the conceptual 
framework differs significantly, in other instances even the mathematical formal- 
ism varies. It is well-known that soon after the inception of quantum mechanics 
de Broglie proposed an alternate formulation based on a guiding wave field, 
an idea rediscovered later by Bohm [i]. Almost at the same time it became 
possible to formulate the new theory in terms of a hydrodynamical model 
and this line of enquiry was later continued and perfected by Takabayasi [3]. In 
both formulations, although in a different context, a prominent place is given 
to what is called quantum potential. The Bohm quantum potential should not 
be considered as another classical potential function acting on physical systems 
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[4]. It is a rich physical concept with broader mereological implications and dis- 
plays with clarity the inherent wholeness and non-separability of an interacting 
many-body system [5] since it does not necessarily diminish with distance and 
it is manifestly non-local. On the other hand its very reliance on a probabilistic 
measure, namely the curvature of the amplitude of the wavefunction [4], seems 
problematic within the scope of a deterministic theory like the de Broglie-Bohm 
one. In the next section we examine some fundamental considerations, introduc- 
ing a non-linear wave mechanical equation and relate its solution with the linear 
Schrödinger equation. In section 3 we examine the hydrodynamic formulation 
and examine the effect of the quantum potential to the osmotic velocity. A time 
symmetrical wave equation with a non-homogeneous term for the probability 
density is also derived that applies to the unmodified dynamical problem. 


2 Foundational issues 


Here we modify the Schrodinger dynamics by introducing an imaginary Bohmian 
potential and then make a polar substitution for the wavefunction. In that way 
we have introduced a quantum potential term in both the Hamilton Jacobi and 
the modified continuity equation as a point of contact and consequently we 
recombine those two expressions in a different manner assuming at the same 
time that the phase of the nonlinear wavefunction is approximately the same 
with that of the linear equation. Our point of departure is a modified nonlinear 
wave equation of the form 
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the Bohm quantum potential discussed earlier. Making a standard polar substi- 
tution of the form ® = Re’S/" and S the multivalued phase of the wavefunction 
gives the continuity equation 
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It follows from (4) that even if the phase does not exhibit any spatial variation 
the probability density may still be time-dependent. It can be reduced to a 
nonlinear equation, transformed to a diffusion equation like (3) and be readily 
solved provided that £ is positive, otherwise the coefficient of diffusion will be 
negative and the problem ill-posed. On the contrary € < 0 describes amplifi- 
cation effects. A diffusion equation requires initial conditions to be imposed in 
order to deal with a well-posed problem, in other words if we use final condi- 
tions for (3) or (4) this engenders the possibility of a blow-up to infinity and a 
singularity in a finite time interval. This means that there is a preferred time 
direction and it is not possible to retrodict solutions in the backward time di- 
rection. If we allow only the first order derivative phase term to remain we 
obtain a diffusion equation with convection that can still be solved exactly. The 
existence of the Bohmian potential in (1) allows us to establish a connection 
between the Hamilton-Jacobi and continuity equations as 
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The right hand of the above expression, which can been seen as a source term for 
the continuity equation, is related to the classical Schrödinger equation [6]. The 
time-independent version of (6) is easily obtained through the relations ce =0, 
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where L the system Lagrangian. Multiplying the last two equations we obtain 
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In that way we have generalized the expression for the quantum hydrodynamic 
propagator given in the de Broglie-Bohm theory [7]. Notice two alternative 
gauge-independent expressions are given. These are formulas relating the wave- 
function ® at two different times and spatial locations, through an exponential 
factor that involves integration over past times. It would be interesting to find 
a relation between ® and W at the same time and location. In order to do this 
we need to make an assumption about the phase. We assume that the phase 
of the unperturbed system does not change significantly when we include the 


non-linear term. It is reasonable to assume that the imaginary term will more 
appreciably affect the evolution equation of the probability density and to a 
lesser extent the time development of the phase since the Hamilton-Jacobi that 
governs the evolution of the phase equation will remain unaltered. Certainly, we 
can add correction terms to this initial phase and use (5) to obtain iteratively 
better approximate expressions. Using (8) and setting € = 0 we find 
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where Y = Aet9S/®. Multiplying both sides with e*%/” and assuming that 
R(ro, to) = A(ro, to) we find 
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Hence we have found a proportionality relation between ® and WV. Notice that 
the quantum potential Qo corresponds to the state Y. This approximate expres- 
sion indicates that the total probability amplitude is rescaled by an exponential 
integral factor with local phase dependency. When the finite duration of the 
time aperture becomes instantaneous so tg + t the exponential becomes unity 
and we recover the linear solution V at equal times. Obviously we recover 
the linear limiting case too when ¢ > 0. Within a many particle description 
the phase modifies a coincidence probability density and features the non-local 
character of the generated dissipation. It must be remembered that an essen- 
tial hypothesis, apart from the inclusion of the imaginary Bohmian potential, 
was the assumption that the phase dependence of the modified wavefunction 
is identical to the linear Schrödinger equation. Only then can we obtain this 
compact relationship between the two wavefunctions which would otherwise be 
unrelated. Since ® is a solution of a nonlinear equation (1) in contradistinction 
to WV it does not satisfy the superposition principle due to the existence of the 
scale factor. Nevertheless, it has the same phase so it may exhibit interfer- 
ence and correlation effects exactly as a coherent superposition of Schrödinger 
equation solutions. A coherent state for example will display the same phase 
properties but an amplitude with a different spatio-temporal modulation. The 
presence of the state associated quantum phase illustrates the kind of indivis- 
ible, essential unity which is characteristic of the whole quantic system, more 
clearly pronounced in the multi-particle case. Normalization concerns dictate 
that ® goes to zero at infinity and it may even involve non-Markovian evolution 
and entail memory effects [8]. Using (12) we can re-express (1) as a non-linear 
Schrodinger equation with an added integral term. This is not a particular 
useful form but it is included here for the sake of completeness 
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Effectively the existence of the imaginary Bohm potential means that the quan- 
tum action becomes complex 
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where a vector potential has been introduced in the usual manner as v = vs = 
<A. The above relation suggests the following form for the non-relativistic 
propagator from place ro, to to r,t 
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This expression clearly implies non-unitary evolution. Again it must be empha- 
sized that it is this phase or action of the system that generates dissipation, so 
it is a kind of self-induced decay and unlike the influence functional or similar 
approaches we do not need to have recourse to a heat bath that affects the sub- 
dynamics of the system under examination. This is of great consequence and 
illustrates the sharp conceptual difference between the way irreversibility ap- 
pears in the present formalism and more common reservoir-system approaches. 
In fact, it may be argued that strictly speaking it is not possible to trace out 
environmental degrees of freedom and divide rigorously a system from its sur- 
roundings, there is no sharp distinction between measured system and environ- 
ment since the quantum phase cannot be appropriated to individual particles 
but points to the irreducible quantum compound and reveals its substantial 
unicity. If the action in the frequency domain is analytic in the upper half 
then its real and imaginary parts are not independent but they are related by 
a Hilbert integral transform 
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the frequency domain representation of the integral scale factor and W(r,w) 
that of (r,t). Furthermore using the Fock-Krylov formula [9] for the survival 
probability of an unstable system we find 


t ae i j 
e$ Sa [E+ (VS—aA)2+V (8.7) ar / g(E)e dE} , (18) 


—oo 





p= 





which in general deviates from exponential decay law even if the density of 
energy states is a Lorentzian function . The decay rate is 
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It is straightforward to generalize (1) for a many-particle system (of equal masses 
for convenience) or even a continuous one like a boson field. It is evident that 
if it is possible to decompose the phase as a sum of functions S = S1 + S2 +.. 
then the multiparticle wavefunction is separable. The rate of change of the total 
probability P = f pdx 
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We see then that the total probability is not conserved but diminishes with time. 
What is important is that (1) entails the possibility of entropy production. Using 
the entropy formula S = — f plnpd?x we obtain 
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which is positive. The positivity of entropy production is a general result fol- 
lowing directly from the modified dynamics and does not rely on the phase 
approximation made earlier. 


3 Quantum hydrodynamics and vorticity 


It was stated in the introduction that quantum mechanics can be expressed 
in terms of quantum hydrodynamic variables, namely the probability density 
and the momentum field, even though the two formulations are not equivalent 
as certain hydrodynamical solutions may not satisfy the Schrodinger equation. 
The inhomogeneous continuity equation for the probability density is analogous 
to the conservation of mass equation in fluid mechanics. It is possible to derive 
an equation analogous to the momentum Navier equation [3] [6] [7]. Employing 
(4) we find the continuity equation for the flux momentum 
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where the quantum force is the sum of the classical and quantum potentials 
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It is possible to write (22) in a tensor form as in [6] 
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where the probability quantum stress tensor is written as 
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It is possible to recover the familiar expression for the tensor equation [6] for the 
limiting case € > 0. Following a development of the Navier-Stokes equations 
found in the aerodynamic sound literature {10} [11] [12] taking the time derivative 
of the continuity equation for the probability density (4) and the divergence of 
(24) and subtracting them we can eliminate the convection term. In that way 
we combine the two in a new exact form as 
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which contains both second and first order time derivatives and in order to find 
its time evolution we need to know initial conditions for p and p, both. The 
second term in (26) contains non-constant coefficients so the equation cannot 
be solved exactly but notice that it resembles the telegrapher’s equation which 
admits solutions expressed by Heaviside step functions. This indicates at least 
qualitatively that the solutions of (26) will propagate not instantaneously but 
with finite speed. As in the aeroacoustics case we subtract from both sides a 
second spatial derivative term multiplied by an arbitrary reference velocity co 
that states the propagation of p and obtain 
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where the driving Lighthill tensor which includes a quantum potential compo- 
nent with highly non-classical properties is expressed as 


(27) 
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For a two-dimensional problem we could write (27) in terms of a stream function 
for the flow velocity v. If we omit the quantum potential term and set £ = 0 in 
(27) we obtain a time symmetric equation which admits solutions in both time 
directions. In contrast, the Schrödinger equation admits solutions only in one 
direction in time and its complex conjugate in the opposite one. In that limit 
(27) is reduced to an inhomogeneous wave equation with a Lighthill source term 
[11] [12] and its solution is 
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where r = |x — y| the distance between a reception point x and the source y. 
Even though this is a formal solution it still allow us to extract some signifi- 
cant conclusions by including both retarded and advanced terms. The existence 
of the retarded term is rather uncontroversial but this cannot be said for the 
advanced term since unavoidably introduces an element of retrocausation |13]. 
A similar issue is well-known in classical radiation theory where the advanced 
terms are excluded on physical grounds, but it is not clear if the same kind of 


argumentation could be employed here. Since the Bohmian quantum potential 
implies nonlocal effects it strongly suggests that we cannot discard advanced so- 
lutions without loss of essential physics. Furthermore, for two or more particles 
the coincidence probability density is connected to the total probability that 
appears in Bell type theorems which is a clear indication that the advanced 
term needs to be included in order to violate the classical correlations. This 
is obviously not a mathematical argument and the matter needs further con- 
sideration but it is plausible to assume that the effect of advanced terms in a 
many-particle context would definitely affect correlations and violate the usual 
assumed conditions for Bell-type theorems [14]. Obviously in the context of a 
classical electromagnetic or acoustic field there is no possibility of violating that 
family of inequalities but this is cerainly not the case for a quantum coincidence 
probability. One of course could escape the conclusion by modifying the ordi- 
nary Schrödinger dynamics. Introducing now the stochastic osmotic velocity 
[7 [15], which is also implied in earlier formulas as in (21) and (25), 
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This is a measurable, phase dependent diffusive velocity guided by the deter- 
ministic flow velocity field. The many particle expression exhibits dependency 
on the remote site locations which is a manifestation of nonlocality. Notice that 
integration over past times can be extended in principle to the initial proba- 
bility distribution which makes a common cause explanation implausible. 
The total velocity is the stochastic one in addition to the standard flux veloc- 
ity. The intrinsic statistical natural of the osmotic velocity field does not allow 
deterministic trajectories and the rule of non-crossing is not applicable. The 
dilatation of the osmotic field is given by 


2 t 
V -u= Vint — = ( 
A Ja 


OV - 

ows) +2(V-v)? +2v-V?2v + vv) dr. (32) 
= 

More interesting is the vorticity of this stochastic velocity field. Taking the 

curl of (31) we confirm that the first term is zero as the gradient of a scalar 

function but the integral term has non-zero contribution and some interesting 

implications. Using the fact that V x v = 0 yields 
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We see then that the quantum velocity has a non-zero component related to 


circulation which could affect in principle the electron current in a conducting 


ring as in the case of mesoscopic Josephson currents. It also follows that other 
quantities related to circulation as the helicity H = fy u-wd*a may be non-zero 
too. The osmotic acceleration is given by 
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Defining the circulation as 
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and applying Stokes theorem we find 
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Hence, the circulation of the osmotic field does change with time and the Kelvin 
circulation theorem is violated. It is important to note that as with the 
entropy production and the tensor equation (27), the non-zero vorticity effect 
is again a general result following from the inclusion of an imaginary quantum 
potential and does not rely on the phase assumption mentioned in section 2. It 
can be immediately verified that the same result holds by substituting the exact 
expression (10) in (30). 


4 Concluding remarks 


We have seen that the addition of an imaginary Bohm potential generates quan- 
tum friction and leads to certain interesting implications. It induces decay ef- 
fects in the amplitude of probability, modifies the exponential decay law, leads 
to entropy production and generates vorticity effects for the osmotic velocity. 
The phase of the wavefunction, and its derivative velocity field, has particular 
importance in the present formalism since it guides the osmotic velocity field. 
Furthermore, still in the hydrodynamic formulation a generalized quantum hy- 
drodynamic propagator has been obtained. In the domain of standard dynam- 
ics, by combining the two hydrodynamical equations a single non-homogeneous 
Lighthill wave equation, manifestly time-symmetric, has been derived which in 
general admits both retarded and advanced solutions. The wider physical im- 
plications of this result suggest that we may need to choose between backward 
causation and exploring the possibility of extending the standard formalism 
beyond the unitary regime. 
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